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1. Introduction 

QCD at a finite baryon chemical potential and/or finite temperature is an important subject of 
study, which is crucial for understanding the early universe, the relativistic heavy ion collisions, 
and the dense matter inside the neutron stars. Although the lattice QCD simulations should play an 
important role for studying the strongly coupled parameter region of this theory, the notorious sign 
problem prevents us from a direct application of the simulation. In principle, a lattice simulation 
can be performed by using the phase-quenched ensemble. The effect of the phase can be taken 
into account by the phase reweighting method. However whether it is practical or not is not clear 
a priori; when the number of the flavors Nf is two, the phase quenched theory is the QCD with the 
isospin chemical potential 1 , whose phase diagram is different from the original theory (for example 
the pion condensation takes place), and hence a severe overlap problem can appear. 

Recently it turned out that the phase quenching and the phase reweighting are actually prac- 
tically useful techniques. The first to emphasized this fact, albeit empirically, are probably Kogut 
and Sinclair [JTJ] _ They pointed out that various model calculations give the same answer for cer- 
tain observables in the full and phase-quenched theories, as long as the pion condensation does 
not take place in the latter. They also pointed out that known lattice data give very similar results. 
Independently, Cohen [Q] and Toublan ^ pointed out the similarity in the large-A^ limit. Recently 
these facts have been understood theoretically unified manner 0,1, |]. Actually there is an exact 
equivalence between the full and phase-quenched QCD at large-iV c , which provides a good approx- 
imation at N c = 3. 2 The equivalence is a version of the large-A^ c orbifold equivalence [^, ||], which 
was discovered through the study of the string theory. (For other interesting applications of the 
orbifold equivalence see e.g. [§].) In this paper, we briefly summarize the equivalence, show the 
lattice data which proves the equivalence can be seen already at N c = 3, and point out the same 
equivalence holds for various effective models in the mean field approximation. 



2. The (partial) equivalence between QCDg and QCD/ in the large-iV c QCD 

2.1 The equivalence between QCD B and QCD/ 

Let us start with introducing the orbifold equivalence. First we choose the discrete symmetry 
P (subgroup of gauge, flavor, or spacetime symmetry) of the parent theory, which is the SO(2iV c ) 
or Sp(2A r c ) theory with the baryon chemical potential (SOg or Spg) in the present case. We then 
throw away all the degrees of freedom not invariant under P. This procedure is called the orbifold 
projection. After the projection, we obtain a new theory called the daughter. We consider two 
different projections, which give QCD with the baryon and isospin chemical potentials (QCDg 
and QCD/) as daughters. The orbifold equivalence states that, in the large-A^ limit, correlation 
functions of operators (A^ , \jf) invariant under P in the parent (called neutral operators) agree 



'As we will see in Sec. 2.2, although this statement is correct for the partition function, there is a difference when 
one considers the physical observables. 

Previously we argued the equivalence is restricted to a class of observables. As we will see, however, the equiv- 
alence holds for any observables. We thank F. Karsch for a valuable critical comment, which made us revisit the issue 
and led to more precise statement. 
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with those of the operators ' (Am , y/" proj ) that consist of projected fields in the daughter: 

{#\ p) e[ p) ■■■) p = (tfWtfW ■ ■ -) d . (2.i) 

Here we take the coupling constants as = g| = g$ p , where the 't Hooft coupling gl v N c is 
kept finite. The field theoretic proof was given in [Jsj] for a class of theories, which can be general- 
ized to various cases. For QCDg, QCD/, SOg and Spg, a couple of evidences of nonperturbative 
equivalence were also provided by the weak-coupling analysis in QCD and QCD-like theories at 
high density limit [||], low-energy effective theories [[!(]], chiral random matrix models [Q] and 
holographic models JTTj] . 

In order to build a projection from SOg to QCD B , we use the Z 4 discrete symmetries of SOg 
generated by J c = —1O2 <8> 1jv c (In is an N x N identity matrix) and CO = e in l 2 e U We require 
the gauge field A s ° ab and the fermion Ya°a t0 t> e invariant under the following Z2 transformation 
embedded in the gauge and U(l)g transformation 

A l°ab = ( J c)aa' A l° a 'b'( J ^ l )b'b, (2-2) 
¥a°a = CO{J c )aa'Wl°a'- (2-3) 

From these projection conditions, it turns out that the daughter is QCD#. The projection symmetry 
breaks down in the BEC/BCS crossover region (diquark condensation region) of SO5, because the 
U(1)b symmetry is broken to Z2 there. 

One can also construct the projection from SOb to QCD/ for even Nf by choosing another Z2 
symmetry [Q. |5|], 

A l°ab = { J c)aa'A S ° cl 'y{ J c V )b'b, (2.4) 
Ya°af = (• / c)oa'Va, a '/'(^r )/'/' (2.5) 

where J { = —ia% <S> generates Z4 subgroup of SU(2) isospin symmetry and the projection 



condition for the gauge field is the same as (|2.2|). In this case, the isospin symmetry used for the 



projection is unbroken everywhere, and so the orbifold equivalence holds including the BEC/BCS 
region of the phase diagram. Therefore, through the equivalence with SOg, we obtain the equiva- 
lence between QCD# and QCD/ outside the BEC/BCS region of the latter; the phase quenching is 
exact for neutral sectors in this region. 

Let us consider the l/N c corrections to QCD^ and QCD/. In the 't Hooft large-A^ limit, 
expectation values of gluonic operators trivially agree because the fermions are not dynamical. 
Now consider finite-A^, say N c = 3 and Nf = 2. Then the largest correction to the 't Hooft limit 
comes from one-fermion-loop planar diagrams, which, as we have seen, do not distinguish jUg 
and /!/. Therefore the difference of expectation values of gluonic operators is at most (Nf/N c ) 
(two-fermion-loop planar diagrams). In particular, the deconfinement temperatures, which are de- 
termined by the Polyakov loop, agree up to corrections of this order. A similar observation was 
made in [||] by a perturbative argument. 

2.2 More equivalence between QCD B and phase-quenched QCD 

It is often said that QCD/ and the phase quenched QCD are the same. However, although this 
statement is correct for the partition function, there is a difference when one considers the physical 
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observables. In QCD/, the propagators of up and down quarks are D 1 (+ju) and D 1 (— fi). On the 
other hand, in the phase-quenched QCD, it is natural to take both D~ l (+fi). (In the terms of the 
lattice QCD simulation, the configuration are generated by using QCD/, while the same operators 
as QCDg are used for the measurement.) Therefore the expectation values of the chiral condensate, 
the baryon density, and the isospin density are calculated as follows: 





QCD B 


QCD/ 


phase-quenched QCD 


chiral condensate 


2<TrD- 1 Gu)) B 


(TrD- 1 (M)+TrD- 1 (-/x)) 7 


2(TrD- 1 (/x))/ 


baryon density 


2(Try°D- 1 ( J u)) B 


(Tr fD~ l (ji) + Try°D- 1 (- J u)) / 


2(Tr7 D- 1 ( J u)) / 


isospin density 





(Tr fD-\n) -TrfD-^-jl))! 






Here ( }g and ( )/ are the expectation values with QCDg and QCD/ ensembles, respectively. We 
can easily see the chiral condensate in QCD/ and the phase-quenched QCD take the same because 
of the charge-conjugation invariance of the QCD/ ensemble. Therefore, the orbifold equivalence 
(the chiral condensate in QCD^ = the chiral condensate in QCD/) tells us it is not affected by the 
phase quenching. For the baryon density, let us remind (Try°D _1 (— ju))/ = — (Try°D _1 (-|-jU))/, 
again because of the charge-conjugation invariance of the QCD/ ensemble. Therefore, the isospin 
density in QCD/ and the baryon density in the phase-quenched QCD take the same value. (Also the 
baryon density in QCD/ becomes zero.) By combining it with the orbifold equivalence (the baryon 
density in QCDg = the isospin density in QCD/), we conclude that the phase quenching does not 
affect the expectation value of the baryon density. The same argument holds for other observables 
too, and the orbifold equivalence leads to the exactness of the phase quenching for any observable 
to 0(N f /N c ). 

Let us provide a heuristic argument which can be generalized to any number of flavors and any 
value of the quark chemical potentials. Let us assume that the overlap problem is not severe, as the 
argument based on the orbifold equivalence shows, and consider why the determinant phase does 
not modify the expectation values. In the large-A^ limit, the quantum fluctuation is suppressed. In 
the phase quenched simulation, with a standard 't Hooft counting, properly normalized operators 
like yy/Nc fluctuates only 0(Nf/N c ). In other words, the histograms of the properly normalized 
quantities have a very sharp peak of the width 0(Nf/N c ). On the other hand, the phase factor is of 
order one 3 , and it cannot change drastically around the peak. Therefore the correction is at most of 
order Nf/N c . 

3. Evidence from lattice simulations at N c = 3 

We have seen that the phase quenching is exact to 0(Nf/N c ). But the standard 't Hooft 
counting does not tell us the expansion coefficients. That motivates us to look at lattice data of 
N c = 3 QCD. In the following we summarize lattice studies which compared QCDg and the phase- 
quenched QCD. 

• In Jl2|], Nakamura et al. studied two-flavor QCDg and phase-quenched QCD by using stag- 
gered fermions. The QCD# is obtained by the phase reweighting. The bare quark mass is 

3 Because the expectation value of the phase factor (e' lmS ) is real, the one-point function (ImS) is zero, and hence 
the leading contribution to the average phase comes from the connected two-point function of the imaginary part of the 
action, {(ImS) 2 ) com ,., which is of order one. We thank Y. Hidaka for a clear explanation on this point. 
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am = 0.05 on a 8 3 x 4 lattice. The chiral condensate and the Polyakov loop are computed for 
a\i = 0.1 and 0.2. They found a perfect agreement between QCD# and the phase-quenched 
QCD within numerical errors. 



• In the right panel of Fig. 1 and the left panel of Fig. 4 of []13[], the free energy at various 
temperatures between 0.5T C and l.ir c are plotted as functions of Q. By putting these plots 
on top of each other, one can see a very nice agreement near the critical temperature and 
<2 < 100. It clearly shows the validity of the phase quenching. It should also be remarked 
that the corrections are still tiny for Nf = 8, a larger number of flavors than Nf = 2 + 1 in the 
real world. 

• Fodor et al. [|l4|] combined the phase reweighting and the density of states methods. In Fig. 4 



of [14] they show the critical couplings at a\i = 0.3 both in the phase quenched and phase 



reweighted cases, which take close values. 

The large-A^ equivalence holds for the imaginary baryon and isospin chemical potentials, 
(XoMrf) = 0'Mim g ,'Aiimg) and {\L u ,\Ld) = O'Mimg, -*'Mimg), without any modification. As a 
result, the chiral condensates {yy)b and take the same value at finite imaginary 

potentials as long as the projection symmetries are unbroken. 



In pig], the chiral critical temperatures T c (jit) in two-flavor QCD were exploited by the ex- 



trapolations from the imaginary chemical potential, by using a fitting ansatz 



They found Q a\ = -0.470(13) for /I/ and a x = -0.522(10) for }X B , which provide a nice 



quantitative agreement already at N c = 3. 

Let us consider the Taylor expansion method, in which the expectation value of an observable 
is expanded in powers of jl/T , 

(0)Bj=ic?-'(>±y. ( 3.2) 

Taylor coefficients c% and which are functions of the temperature T, can be determined 
by the simulation at /I = 0. The large-A^ equivalence tells that the coefficients agree in the 
large-A^ limit: lim^^ooC^ = lim^-j-ooC^. 

In [15 1, the coefficient cf and c\ for the chiral condensate and the pressure of the quark- 



gluon gass have been calculated in two-flavor QCD. Although the difference between cf 
and c ! 2 are not very small for T < T c in the chiral symmetry broken (and confinement) phase, 
they agree exceptionally well for T >T C . The origin of the difference for T <T C may come 
from the contributions of thermally excited mesons which large-A^ is suppressed at large- 
N c . On the other hand, for T > T c , fundamental degrees of freedom are deconfined quarks 



4 For odd n, cf and c ! n vanish, and the first nontrivial /x -dependences appear in cf and c\. Although cf (n > 4) have 
been calculated, c/ n (n > 4) have not been calculated in [|l5|. (Note that, for n > 4, they use the same symbol (/ n for 
another quantity.) 
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and gluons rather than baryons or mesons, where the difference between QCDg and QCD/ 
becomes much smaller and the large-A^ equivalence is very well satisfied even at N c = 3. 

4. The equivalence in the effective models 

As an example of the effective models, we consider the Nambu-Jona-Lasinio (NJL) model. In 
order to simplify the discussion, we consider the chiral limit. The starting point is the Lagrangian 
with the \J(N C ) color current interaction with Nf flavors, 

j^njl = ¥f (fd,+ m 4 ) Yf + jftSrf!! > (4-D 

where Jjj^ = WfY^^vYf an d Tfj w& tne U(-/V c ) color generators and summation is taken over 
repeated indices. The coupling constant G is taken to be of order A^ . One rewrites it keeping only 
the interactions in the scalar and pseudoscalar channels after Fierz transformations: 

^njl = +/!//) Y + ^int, 

■%nt = jj- [(VfVf)(VfVf) + (Wf^Wf){Wf^Wf)] • (4-2) 



In the Lagrangians ( |4.l| ) and ( |4.2[ ), the invariance under U(A^ C ) gauge symmetry and U(Afy)z. x 
U(N/)r flavor symmetry are manifest. Here we ignore the effect of instantons or the U(1)a 
anomaly which explicitly breaks the U(1)a symmetry, because it is subleading in l/N c . (From 
the viewpoint of the orbifold equivalence, there is no reason for the exactness of the phase quench- 
ing at the level of mean-field approximation (MFA) if we take into account the 1 /A^-suppressed 
instanton effects. However, even if we incorporate them, the phase quenching for the chiral con- 
densate turns out to be exact within the NJL model [[J.) For SO^A^.) theory, we can construct 
the corresponding NJL model in the same manner, by using the S0(2N C ) current. The proof of the 
large-A^ orbifold equivalence applies to the NJL model, by starting with the NJL model for SOg 
and by using similar projection conditions as the previous section 

When one considers the large-A^ limit, one should set up the correct 1 /./Vc -counting scheme 
which reproduce the correct 1 /A^-scaling in the large-A^ QCD. The quark y has N c colors so that 
a closed color loop gives a factor of N c . The coupling constant of the four-fermi interaction should 
be taken as 0(N ( 7 l ), and furthermore, the form of possible four-fermi interactions are restricted; in 
other words only the interactions which have origins in QCD are allowed. In this setup, the right 
1 /./V c -counting follows and we can use the same proof of the orbifold equivalence as the large-A^ 
QCD [§]. 

Now let us see the relationship between the large-A^ limit and the MFA. We first perform the 
Hubbard-Stratonovich transformation by introducing auxiliary fields corresponding to the fermion 
bilinears, a = and % a = 1//7/5 % a y, and then integrate out fermions to obtain the partition funcion 

Z = e- W = [ dodne- I{ - a ^ . (4.3) 



Here I(a, n) is the bosinized effective action 



r(a,K)=N c 



-TrlogD + ±Jd 4 X (o 2 + n 2 a ] 



(4.4) 
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with D = y^d^ + 2oa + 271a- It describes mesons a and %. Because of an overall factor N c in the 
action, the 1 /N c -expansion is equivalent to the expansion with respect to meson loops. The leading 
order conesponds to the saddle-point approximation, or equivalently the conventional MFA where 
the auxiliary fields are replaced by the expectation value (i.e., the mean-field). In order to go beyond 
the MFA, we have to take into account meson-loops order by order. 

Similar arguments hold also for various other theories, such as linear sigma model Polyakov- 
Nambu-Jona-Lasinio model [0], Polyakov-quark-meson model [0|, chiral random matrix model 



[g], Sakai-Sugimoto model []6| and D3/D7 model [jll|]. 

5. Conclusion and Outlook 

We have seen the phase quenching is exact to 0(Nf/N c ), outside the pion condensation of 
the phase quenched theory. In other words, the effect of the phase is 1 /N c -suppressed, and hence 
the reweighting method works without being threatened by the overlap problem. Previous lattice 
studies confirm the effect of the phase is small already at N c = 3. We have also shown the exactness 
of the phase quenching in effective models, which had been realized by explicit calculations and 
used to justify the reweighting method, can be understood in a unified manner, from the point of 
view of the large-A^ equivalence. Now the phase quench and phase reweighting methods have a 
theoretical justification, and hence it is important to study the QCD phase diagram by using them. 
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